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ABSTRACT 

The nonlinear evolution of rclativistic magnetic reconnection in sheared magnetic configuration 
(with a guide field) is investigated by using two-dimensional rclativistic two-fluid simulations. Rcl- 
ativistic guide field reconnection features the charge separation and the guide field compression in 
and around the outflow channel. As the guide field increases, the composition of the outgoing energy 
changes from enthalpy-dominated to Poynting-dominated. The inertial effects of the two-fluid model 
play an important role to sustain magnetic reconnection. Implications for the single-fluid magneto- 
hydrodynamic approach and the physics models of relativistic reconnection are briefly addressed. 
Subject headings: magnetic fields — MHD — plasmas — relativity 



1. INTRODUCTION 

The role of magnetic fields has been recognized in var- 
ious contexts of relativistic astrophysics: pulsars, mag- 
netars, gamma-ray bursts (GRBs), active galactic nuclei 
(AGNs), and black holes. Owning to its fundamental 
nature, relativistic magnetic reconnection in electron- 
positron pair plasmas (or in pairs and baryons) has 
drawn attentions in these sites as well; however, its mech- 
anisms are poorly underst ood. Despite several theoreti- 
cal attempts over d ecades (Blackma n fe Field|1994 Lyu 



tikov & Uzdensky 2003 Lyubarsky 2005}, the most im- 
portant properties of relativistic magnetic reconnection 
such as the energy conversion rate and the released en- 
ergy composition are still under debate, especially in the 
magnetically dominated limit. 

Relativistic magnetohydrodynamic (RMHD) codes are 
very desirable to study the theories of relativistic mag- 
netic reconnection and to investigate astrophysical prob- 
lems which involve magnetic reconnection. However, 
surprisingly, reconnectio n or non-ideal RMHD problems 
were left untouched until |Watanabe fe Yokoyama ( 2006 1 
carried out resistive RMHD simulations. TEey exhib- 
ited a Petschck-like structure with an Alfvenic outflow 
and showed that relativistic reconnection is faster than 
the nonrelativistic counterpart. Several groups have de- 
veloped resistive RMHD codes (|Komissarov|2007| |Palen- 



zuela et al.|2009||Dumbser fe Zan otti 2009}, which can be 
applied to the r econn ection problems as well. Recently, 



Zenitani et al. (2009) employed a relativistic two- fluid 
model to simulate magnetic reconnection. They obtained 
Petschek-type steady reconnection in a large system and 
found that the reconnection speed becomes faster and 
faster in magnetically dominated regimes. 

In addition to the basic antiparallcl configuration, re- 
connection with a current-aligned magnetic field (the 
"guide field" ) is likely in magnetic shear and in celestial 
flare situations. It is well known that guide field recon- 
nection behaves differently from its anti-parallel counter- 
part. In the relativistic regime, it is expected to change 
the energy compo sition in the outflow in RMHD scales 



( |Lyubarsky||2005| ). However, the RMHD/fluid-scale be- 
havior of relativistic guide field reconnection has not yet 
been explored by simulations. 



In kinetic scales, it has been found that reconnection 
is a powerful particle acceler ator in both antiparallel 
and g uide field configurations |Zenitani fe Hoshino|2001| 
20071 20081 IJaroschek et al.|20(M]r^ mce the plasma non- 
thernal energy is comparable with or exceeds the thermal 
part, we expect that kinetic physics significantly affects 
the global dynamics. Comparison with RMHD/fluid 
models is useful to understand the role of kinetic physics. 

The purpose of th is paper is to advance our two-fluid 
reconnection work ( Zenitani et al.|[2009[ ) another step 
forward. In the earlier work, we assumed symmetric 
motions of the two species, which enforced the charge 
neutrality in the system. In this work, we solve two fluid 
motions independently. This allows us to explore broader 
ranges of physical targets, which involve charge separa- 
tion. Using relativistic full two-fluid simulations, we will 
investigate a nonlinear development of a relativistic re- 
connection system with a guide field. 

2. NUMERICAL SETUP 

We use the following set of relativistic fluid equations 
and Maxwell equations. The subscript s denotes species 
('p' for positrons and 'e' for electrons). Similar equations 
for electrons are considered too: 



d_ 
Ft 



l p n p = -V • (ripUp), 
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TABLE 1 
List of Simulation Runs 



In these equations, 7 S is the Lorentz factor, n s is the 
proper density, u s — j s v s is the fluid 4- velocity, w s = 
n s mc 2 + [r/(r — l)]p s is the enthalpy with the specific 
heat ratio T = 4/3, Sij is the Kronecker delta, p s is 
the proper pressure, q p — —q e is the positron/electron 
charge, j = q p n p u p + q e n e u e is the electric current, and 
Pc = lp n p% + 7e^e9e is the charge density. In the mo- 
mentum and the energy equations, inter-species friction 
terms with the coefficient Tf r are added. They are re- 
vised from our previou s work. The new form is s imilar 
to a covariant form in Gurovich & Solov'ev ( 1986 ) . The 
variables i/> and <j> are virtu al potentials for dive rgence 



cleaning (Munz et al. 2000 Dcdncr et al. 2002 Komis 



sarov|2007| , whose propagation speed and the decay rate 



are c and k, respectively. 

We study two-dimensional system evolutions in the 
x-z plane. We use the Harris-like model as an initial 
configuration: the magnetic field, the (proper) density, 
and the pressure are B(z) = £? tanh(z/L) x + B G y, 
n s (z) — no cosh~ 2 (z/L) + rij n , and p s (z) — n s (z)mc 2 , 
respectively. Here, L is the current sheet thickness, Bq 
is the guide field amplitude, and rii n = O.lrio is the back- 
ground proper density. The initial out-of-plane current 
is sustained by the fluid drift of ~ ±0.1c. This restricts 
the electron inertial length (~0.1L) and the Larmor ra- 
dius (<0.05L) to sufficiently small values. We set the 
reconnection point at the origin (x, z) ~ (0, 0) in the 
system domain of [0, 120L] x [-60L, 60L] or [0, 240L] x 
[— 60L, 60 L]. In order to reduce the computational cost, 
we assume that the system is point-symmetric around the 
reconnection point: f(x, z) — f(—x, —z) or —f(—x, —z) 
for all times. The boundary conditions at x — are set 
accordingly. At the other three boundaries, we basically 
consider the Neumann boundary conditions (df/dn = 
0), and the normal components of the fields (B n and 
E n ) are adjusted to satisfy V • B — and V • E = 47rp c . 
We employed the spatially localized resistivity, which is 
controlled by the coefficient rt r . The effective Reynolds 
number of the frictional resistivity is S = 30 near the 
reconnection point and S = 3000 in the background re- 
gion. We normalize timescales by the light transit time 
t c = L/c. The decay time scale for divergence cleaning 
is set to = 0.2t c . 

We studied the guide field effect by changing Bq as 
shown in Table [ll We consider the magnetization pa- 
rameters a = B^[4tt(2^ 2 w)], the ratio of the Poynting 
flux to the plasma energy flux. In the inflow region, 
<J%n ~ cr x ,in + Vy,in> where a x = B 2 . / [4tt (2j 2 w)] and 
a y = B 2 /[4tt(2j 2 w)] are contributions from B x and B y . 
The parameter a x ,in is set to 4 and the relevant Alfven 

velocity is c A ,in = [<J x ,in/Q + cr ^,m)] 1/2c = 0.894c. The 
other parameter a y ^ n depends on Bq. 

The system evolution is solved by the modified Lax- 
Wendroff scheme with a small artificial viscosity, whose 
viscous coefficient depends on the local gradient of the 
fluid 4- velocity. We directly s olve quartic equations to re- 
cover the primitive variables ( Zenitani et al.| [2009). The 



Run 


Domain Size 


Grid Points 


B G /Bo 




■R 


1 


120 X 120 


2400 X 2400 


0.0 


0.0 


0.142 


2 


120 x 120 


2400 X 2400 


0.25 


0.25 


0.126 


3a 


120 x 120 


2400 x 2400 


0.5 


1.0 


0.102 


3b 


240 x 120 


4800 x 2400 


0.5 


1.0 


0.102 


4 


120 x 120 


2400 x 2400 


1.0 


1.0 


0.074 


5 


120 x 120 


2400 x 2400 


1.5 


9.0 


0.055 


Note 


— The initial guide field 


Bq , the 


ratio of 


the its 



Poynting flux to the plasma energy flux cr y 
reconnection rates TZ. 



and typical steady 



results are checked by changing the system size, the left 
boundary conditions (with or without the point symme- 
try), the time and spatial resolutions, and parameters for 
the numerical stability. 

3. SIMULATION RESULTS 

3.1. Overview 

Snapshots of runs 3b with a guide field Bq/B q = 0.5 
are presented in Figure [ll Shown in Figure [lb is an 
outflow structure of the plasma mean number flow, 



(v = 



0p Uj p 



l p n p + 7 e n e 



(8) 



We see that a fast outflow jet travels to the magnetic 
island (plasmoid) in the downstream region. Since there 
are ambient current sheet plasmas, the plasmoid exhibits 
a crab claw structure in the downstream, similar to large- 
scale MHD si mulations of nonrelativ istic magnetic recon- 
nection (e.g., Abe & Hoshino (2001 1). Inside the outflow 
channel, the flow speed {v x } is 0.8 ~ 0.85c. The max- 
imum 4- velocity of each species is u xs ~ 2.1c. These 
speeds are slightly slower than the antiparallel counter- 
part. The backward flows around x ~ 70L (blue re- 
gions; (v x ) < 0) are reverse flow from the plasmoid. 
Since the downstream plasmoid expands, the surround- 
ing magnetic fields are compressed and then its magnetic 
pressure expels the plasma along the field lines. Sharp 
boundaries at x ~ 60L are the fronts of such backward 
flows. 

Shown in Figure [lj> is the out-of-plane magnetic field 
By /Bo. We see that the guide field is compressed at 
the edge of the plasmoid and in the narrow reconnec- 
tion outflow region, because the incoming reconnection 
flows transport and pile up the B y magnetic flux from 
the upstream region. Such Sy-compression changes the 
composition of outgoing energy flux in relativistic mag- 
netic reconnection. 

Figure [lb shows the charge separation, (j P n p — 
"f e n e ) / {"f P n p + 7e^e)- The charge density p c looks simi- 
lar. We see that positrons dominate in the upper of the 
outflow region, and electrons in the lower. The charge 
separation generates the electric fields in the x-z plane: 
E x < in the upper half of the inflow region, E x > 
in the lower half of the inflow region, and E z < in the 
rightward outflow region. The E x B condition with the 
guide field B y is consistent with the reconnection flow 
pattern. It is impressive to see that the non-neutral lay- 
ers extend over the large spatial domain and that the 
charge separation is strong ^0.5. Such a non-neutral 
structure is qualitatively consistent with positron-rich 
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Fig. 1. — (Color online) Snapshots of run 3b at t = 200r c in the x-z two-dimensional plane. Contour lines show the in-plane magnetic 
fields, (a) The average plasma outflow (v x ) = (npU XiP + n e n a ; I e)/(7pip + 7e^e)> (b) the out-of-plane magnetic field B y /Bo, and (c) the 
charge separation (l v n v — 7 e rt e)/(7 P n p +7 e n e ). (d) The out-of-plane magnetic field B y /Bo in run 3b at t = 400t c . We discuss properties 
along the white line in Section |3.4| 



or electron-rich regions around the reconnecting diffu- 
sion r egion in previous kinetic work ( jZenitani fc Hoshino j 
2008 1. 1 Stripes around x < 60L are transient oscillations 
invoked by the reverse flow. As the system evolves, the 
plasmoid and the reverse flows move further away. We 
also see strong perturbations around x ~ 120L, where 
the plasmoid hits the ambient plasmas. 

Later, the system approaches to a quasi-steady stage. 
Shown in Figure [Tlf is a late-time structure of B y at 
t = 400r c in run ah. We see that the outflow struc- 
ture in Figure [TJ> straightforwardly extends toward the 
+x direction. The backward flow fronts are outside the 
presented domain at this time. The in-plane magnetic 
fields show a narrow Petschek-type structure. Interest- 
ingly, the By-structure is weakly bifurcated, and then the 
two peaks separate the outflow channel into three layers. 
We examine the outflow channel structure in Section 13.41 
Around x ~ 115L, the bifurcated peaks exhibit a very 
weak wave structure, probably due to a velocity-shear 
driven instability. This introduces a minor oscillation in- 
side the central channel in the further downstream; how- 
ever, it does not change the global outflow structure. We 
also confirmed that the global outflow structure does not 
depend on the domain size by comparing runs 3a and 
3b. Therefore, the outflow boundary condition does not 
change the physics. 

We believe that the large distance to the upstream 
boundary (z = 60L) is virtually eliminates boundary ef- 
fects on the simulation. The plasmoid and reconnected 
field lines are well confined in \z\ < 10-20L in panels in 
Figure[l] In the steady stage of t = 400r c , the slope angle 



of the field lines is only ~ 0.1 around the outflow region, 
and it becomes even smaller further upstream. For ex- 
ample, in the vicinity of the top boundary, the field line 
at (0,+56.2L) is connected to (120L, +60L). 

3.2. Reconnection electric field 

The reconnection electric field E y at the X-point 
presents the transfer rate of the antiparallel magnetic 
fields from the upstream region. Since the surrounding 
conditions change in time, the normalized rate or "recon- 
nection rate" gives a good measure of the system evolu- 
tion. Here we define the rate as 



K 



cE„ 



CA-. 



\B~ 



(9) 



1 Charge sig ns are opposite from Figure 2(c) in [Zenitani 
|Hoshino| ( |2008 I , because the guide field is set oppositely! 



where the subscript in' denotes the upstream properties 
at x = 0, z = 20L. The time evolution of the rates is 
shown in Figure [2)3, and we see that the (normalized) 
reconnection becomes constant after the initial ramp up 
for t > 50t c . For comparison, the rate for run 3a is 
presented, too. Runs 3a and 3b are in excellent agree- 
ment except for around t ~ 400r c (indicated by an arrow) 
when the boundary effect temporally comes back to the 
reconnection point in the smaller run 3a. In Table [T] we 
present the peak values of the quasi-steady reconnection 
rates for all other runs. 

So, what sustains such quasi-steady evolution of mag- 
netic reconnection? In other words, what is responsible 
for the reconnection electric field E y around the X-point 
in our simulations? We study the composition of E y near 
the X-point using the Ohm's law. Combining the mo- 
mentum equations (Equation [2]) of the two species, we 
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(a) Reconnection rates (b) Composition of Ey 




Light Transit Time (L/c) Inflow Position (z/L) 

Fig. 2. — (a) Time evolution of reconnection rates: normalized rates 7Z for runs 3a and 3b, and a raw reconnection rate for run 3b 
Ey/Bo. (b) Composition of the reconnection electric field E y along the inflow line (x = 0) in run 3b at t = 200r c : the convection term, 
the advective inertia, the frictional resistivity, and the viscous inertia, as presented in Equation [TT] These values are normalized by Bq. 



obtain the following relation: 



x B 



1 



c (l P n p + r y e n e )q p 

m^ p n p (v p ■ V)hpU p - m~f e n e (v e ■ V)h e u e 

+ (Vp p - Vp e ) + 2Tf r n p n e (u p - u e ) , (10) 



where h s — w s /(n s mc 2 ) is the specific enthalpy. Note 
that we drop the time derivatives considering the quasi- 
steady condition (d t ~ 0). In addition, along the inflow 
line (x = 0), we find n p = n e , j p = y e , v z , p = v z>e = (v z ), 
hpU V:P = —h e u V} e, and B z = 0, and so the equation can 
further be simplified. We can also approximate that the 
artificial viscosity works for h s u s in a form of —v z d zz 
with an effective viscous coefficient v z . The ^-derivative 
terms (d x and d xx ) are negligible in this case. As a result, 
the y-component of Equation [10] along the inflow line 
yields 



E,, 



{v z )B x m(v z ) dhpU y .p r) eff . 
c q p dz 7 P q p 



mv z d hpUy^p 



dz 2 



(11) 



where T) e ff = (jfrjq^) is the effective resistivity by the 
inter-species friction. Terms in the right-hand side rep- 
resent the field convection, the advective inertia of the 
fluid, the frictional resistivity, and the viscous inertia, 
respectively. 

Based on Equation |11| we study the composition of 
the reconnection electric field E y in Figure [2fc. Here, v z 
is represented by the typical value aroundthe neutral 
point. Outside the central diffusion region, the convec- 
tion term mostly explains the electric field. The frictional 
resistivity term works inside the diffusion region where 
the current is strong. However, it explains only a quarter 
of E y . We find that the local Lorentz factor 7 P ~1.6 par- 



tially suppresses this term in Equation 11 Instead, the 
other two incrtial effects are important. The advective 
term is the biggest contributor. In addition to the accel- 
eration in the ±y directions (u VlS ), the specific enthalpy 
h s ~ 1 + Ap s /(n s mc 2 ) increases because dissipation pro- 
cess heats incoming plasmas. In a very vicinity of the 
neutral plane (z = 0), the viscous inertial effect replaces 



the advection. Although we introduced it for the nu- 
merical stability, we believe it reasonably represent the 
physics because the kinetic effect plays a qu asi- viscous 
role near the neutral plane fHesse et al.|2004 l. 

We carried out a similar analysis in the antiparallel 
case (run 1), too. The frictional resistivity plays a bigger 
role, but it explains only one-third of the electric field. 
The local Lorentz factor 7 p ~1.6 similarly suppresses the 
frictional resistivity. Regardless of the presence of the 
guide field, the inertial effects play a role to sustain mag- 
netic reconnection in the relativistic two-fluid system. 

3.3. Energy budget 

Next, we investigate how the guide field changes the 
energy conversion process near the reconnection region 
in the quasi-steady stage. We consider the composition 
of the energy flux F in the following way: 



cE x B x - 

s—p,e 

cE x B y y cE x (B x x + B z z) | ^ ^lsP& 



(12) 



47T 

„2 \ " 



4tt 



E 



r-i 



+mc 2^(7s - l)n s u s + 



(13) 



In Equation 13 we decompose the Poynting flux into the 
contribution by the guide field (B y ) and the one by the 
reconnecting magnetic fields {B x and B z ). The plasma 
energy flux is decomposed into three parts. The pressure 
term transports the gas internal energy and it also con- 
tains the work by the gas. Since it recovers the classical 
enthalpy flux in the nonrelativistic limit (— > ^p s v s 
with T = 5/3), we call this term "enthalpy flux" in this 
work. Note that the rest-mass contributions are consid- 
ered separately from this enthalpy. The last two stand 



for the bulk kinetic energy (— > ^ s 



,v 3 ) and the 



matter flow (- 



mc 



^2 s n s v s ) 7 respectively. 



We consider the energy budget around the reconnect- 
ing square region of \z\ < 10L, \x\ < 40L. The top panel 
in Figure [3] presents the composition of the incoming en- 
ergy flux (— F ■ z) per cross section at the inflow bound- 
ary (z — 10L) as a function of the guide field. They are 
normalized by the typical Poynting flux c(B 2 /4tt) at the 
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upstream border x — 0, z = 10L. The bottom panel in 
Figure [3] presents the outgoing energy flux (— F ■ x) at 
the outnow boundary (x = 40L), evaluated by the same 
unit as the upstream values. They are evaluated at the 
specific time steps in different runs, but they are care- 
fully selected so that the energy budget becomes nearly 
steady in the rectangle region. We see that both two 
fluxes are well balanced. Although it may be difficult to 
recognize in the figure, the matter fluxes are also well 
balanced. 

For comparison, rescaled values of reconnection rates 

t'able [TJ are overplotted on the top panel in Figure 
It is reasonable to see that they are proportional to 
e incoming energy flux except for _Bj,-Poynting flux. 
The reconnection rate is relevant to the inflow speed vi n , 
which transports the upstream antiparallel magnetic flux 
B x into the reconnected field B z . In these cases, under 
the similar upstream conditions, the i? x -Poyting flux de- 
creases as and the accompanying plasma flux also 
decreases accordingly. We also see that the reconnection 
rate decreases but it weakly depends on the guide field 
amplitude Bq/Bo. This trend is cons istent with a lot of 
previous nonrelativistic surveys (e.g., |Huba| ( |2005] |). 

We see that reconnection well dissipates the upstream 
antiparallel magnetic energy. When the guide field is 
weak, the Poynting energy is mostly converted into the 
plasma energy. Importantly, the energy is converted to 
the enthalpy flux rather than the bulk kinetic energy. 
The plasma pressure becomes strong in the outflow re- 
gion in order to balance the strong upstream magnetic 
pressure. In magnetically dominated cases, plasma tem- 
perature becomes relativistic, and then such relativis- 
tic pressure substantially enhances the enthalpy flow. 
Therefore, the enthalpy flux dominates the bulk kinetic 
energy. 

As the guide field increases, the relevant Poynting 
flux becomes a major component of the incoming en- 
ergy flux. Simultaneously, the outgoing energy flux is 
dominated by the Poynting flux of the guide field. In 
the case of Bq/Bq = 0.5 (run 3), as presented in Sec- 
tion 3.1, the guide field is compressed by a factor of ^3. 
Consequently, the outgoing guide field Poynting flux in- 
creases and then it exceeds the enthalpy flux. When 
the guide field is stronger, the outflow is dominated by 
the Poynting flux by B y . In the outflow region, the ra- 
tio of the By-Poynting flux to the plasma energy flux is 
& v ,out — 0, 0.32, 1.21,4.29, and 8.99, respectively. It is in- 
teresting to see that these ratios are similar to the initial 
upstream a y .i n parameter (Table [l]). 

Under the single-fluid ideal RA1HD approximation, we 
have the following relations from the polytropic law, the 
continuity, and the out-of-plane flux conservation: 



d_ 
dt 



0, ±(^L 
dt\j'n' 



0. 



(14) 



where the prime sign (') denotes the single-fluid proper- 
ties. They immediately suggest that a y increases through 
the weak compression by magnetic reconnection, 



din <j y 
dt 



dt 



(2-r) 



dhir 
dt 



(15) 



In the case of run 3b, the typical compressional ratio 
of ~ 2 suggests that a y increases by a factor of ~ 1.6. 



However, since our two-fluid system contains non-ideal 
effects, j^(p' /n' v ) > tends to reduce o~ y . In the stronger 
guide field cases, it is known that plasmas behave incom- 
pressibly and so the system is more likely to preserve 

0~y,in ^ O~y,out- 

Incoming energy 




Bq/Bq 



Outgoing energy 



Poynting flux (By) 
Enthalpy flux 
Matter flow 



■ Poynting flux (Bx, Bz) 

■ Bulk kinetic energy 

-•- Reconnection rate (rescaled) 



Fig. 3. — (Color online) The incoming and outgoing energy fluxes 
around the reconnection region (|z| < 10L, \x\ < 40L). The guide 
field Poynting flux, the rest part of Poynting flux, the plasma en- 
thalpy flux, the bulk kinetic ene rgy , and the matter flow are pre- 
sented, as presented in Equation |13| The reconnection rates 1Z are 
rescaled from Table [T] The relevant a y parameter is indicated by 
small numbers on the top/bottom of the bars. 



3.4. Outflow channel 

We look at the structure of the outflow region. Shown 
in Figure [4] are one-dimensional profiles across the out- 
flow channel, at x = 90L at t = 400r c . The cut line is in- 
dicated by a white line in Figure [TJf. From these profiles, 
we see that the outflow region consists of three character- 
istic layers, separated by two peaks of B y : (1) a positron- 
rich boundary layer on the upper side (2L < z < 4L), 
(2) a central channel with high plasma pressure, and (3) 
an electron-rich boundary layer on the lower side. 

The central channel contains a minor oscillating struc- 
ture in plasma pressure (Figure [4]a) . Inside the central 
channel, the plasma temperature is hot and the outgoing 
fluid velocities are roughly similar v S:X ~ 0.8c across the 
channel. Since the Lorentz force q s (v SjX B y ) works differ- 
ently in the presence of the guide field B y , the positron 
and electron outflow channels are oppositely displaced 
to the indirections. The structures of the density, the 
pressure, the outflow speed v x , and the inflow speed v z 
(magnified 10 times larger in Figure 4]c) are all consistent 
with the displacement of the flow channels. The central 



ZENITANI, HESSE, & KLIMAS 



channel becomes wider in the further downstream and we 
confirm that the system keeps the three-layer structure 
at least x ~ 180L. 

We see that currents in the boundary layers are respon- 
sible for the global magnetic field topology. For example, 
j x explains the compression of the guide field B y , and a 
"bifurcated" j^-structure is consistent with the Petschck- 
type structure of in-plane magnetic fields (Figure [T]f). 
We also see minor reverse currents near the center, Dut 
they are less important. Note that relative motion be- 
tween the two species is significant in the boundary lay- 
ers (Figure [4)2) , because plasmas sustain the currents in 
low-dense regions. Interestingly, in the upper positron- 
rich boundary layer, we see a fast electron flow in the 
y-direction and indeed electrons are major contributors 
to j y there. This is because the Lorentz force j x B z mod- 
ulates the local bulk outflow to the — y-direction, and so 
the electron flow looks overemphasized in our simulation 
frame. 

In addition, in order to see the difference between our 
two-fluid model and the single-fluid RMHD model, we 
evaluate an equivalent single-fluid properties in the fol- 
lowing way. We translate the conserved properties with 
the single-fluid properties (with the prime sign), 

jpn p + j e n e = j'n', (16) 

jpW p u p + "f e w e u e = ^'w'u'. (17) 

7pWp -p p + 7gU> e ~p e = j' 2 w' - p' , (18) 

and then we calculate the sing le-fluid primitive var iables 
by solving a quartic equation ( Zenitani et al.|2 009). The 
single-fluid pressure is assumed to be isotropic. Com- 
bining two fluid pressure into a single-fluid pressure does 
not lead to an isotropic pressure unless the relative ve- 
locity vanishes. The obtained pressure p' is overplotted 
in Figure [4]a. We see that p' is in excellent agreement 
with the total two-fluid pressure (p p + p e ) in the central 
channel. On the other hand, the discrepancy is signifi- 
cant p' I (p p + p e )^lA in the two boundary layers and so 
it may widen the outflow channel in single-fluid simula- 
tions. It is reasonable that the discrepancy is observed 
where the relative motion between the species is signifi- 
cant (Figure |4]c) . 

In the boundary layers the local frozen-in condition 
also breaks down. We similarly analyze the composi- 
tion of t he rcconnection electric field E y based on Equa- 
tion 10 Figure Hlff presents the profile of E y , along with 
the held convection and the z-convection of fluid inertia 
terms, 
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Fig. 4. — (Color online) Structure of the outflow channel at 
x = 90L at t = 400t c . (a) Proper pressure p p , p e , positron 
density n p , n e , total plasma pressure p p + p e , and the relevant 
single-fluid pressure p' are presented. (6) The out-of-plane mag- 
netic field By, and the electric currents j x and j y . (c) Fluid ve- 



locities, v x 



»ij»"iei%p:%ei%pi an d t) Je . The ^-components are 
magnified 10 times larger, (d) Composition of the electric field, E y , 
f— lv ) X B/c)„. and the z-convection of the inertial terms (Equa- 
tions T19] and |20|l. The inertial terms are magnified 10 times larger. 



m-f p np 



Qpilprip + 7e^e) 



OhpUy p 



0h e Uy e 



dz 



(19) 
(20) 



In Figure [4jf, the inertia terms are magnified 10 times 
larger to emphasize them. We found that the non-ideal 
contribution (the difference between E y and the (v) x B 
term) is mostly explained by these terms. The positron 
inertial term appears in the lower electron-rich bound- 
ary layer, and electron term in the lower electron-rich 
layer. For example, in the positron-rich boundary region, 
both the z-component of the positron velocity and the y- 
momentum of the positrons are much smaller than that of 



the electrons. As a result, the electron inertial term dom- 
inates even though the positron density is larger. Around 
z ~ ±3L, the inertial contributions from both two fluids 
sustain 12% of E y . Because of the large Reynolds num- 
ber S ~ 3000, the frictional resistivity effect is negligible 
and we confirmed that the viscosity is negligible, too. 
Compared with the reconnection site where the nonideal 
terms are 100% responsible for the electric field, the non- 
ideal effects are smaller, but they still play a role here. 

Shown in Figure [5] is the x-dependence of the ratio 
of the single-fluid pressure (p') to the total fluid pres- 
sure (p p + p e ), which is a good measure of the two-fluid 
effects. As the central outflow channel expands down- 
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stream, the two-fluid effects are always localized in the 
relevant boundary regions. Their amplitude and spatial 
width are similar until the backward flows around the 
plasmoid hit the boundary layers. 
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Fig. 5. — One-dimensional profiles of p'/{p P +p e ), the ratio of 
the single-fluid pressure to the total two-fluid pressure. They are 
taken from x = 30, 60, 90, and 120L at t = 400t c in run 3b. 



4. DISCUSSION AND SUMMARY 

Regarding basic models of antiparallel reconnections 
there have been two contro versial opinions on the rel- 



ativistic reconnection rate. |Blackman fc Field ( 1994 1 
argued that relativistic reconnection can involve rela- 
tivistic inflow and the fast consumption of the upstream 
magnetic energy, because Lorentz contraction of the 
plasma outflow enables large r energy output p er cross 
section. On the other hand, Lyubarsky (2005J pointed 
out that relativistic pressure increases an effective in- 
ertia w, and that the shock balance condition allows a 
narrower Petschek outflow. Therefore, he argued that 
relativistic reconnection will not be an efficient energy 
converter, due to the slow outflow and the narrower en- 
ergy output channel. Numerically, Zenitani et al. (2009) 
demonstrated that reconnection tends to be faster and 
faster in magnetically dominated regimes, although the 
obtained Petschek -type structure is best described by the 



Lyubarsky ( 2005 1 model. The unsolved problem is what 
makes reconnection faster, or what balances faster energy 
input? 

In our energy balance analysis, we found that the out- 
going plasma energy is mostly carried by the enthalpy 
flux or the internal pressure flux in the antiparallel case. 
We think this is an important reason why relativistic 
magnetic reconnection is faster. The plasma pressure 
in the outflow region needs to be relativistic in order to 
balan ce a strong upstream magnetic pressure ( Lyubarsky 
2005). As the outflow plasma temperature becomes rel- 
ativistically hot, we see that the enthalpy flux (pres- 
sure part) vastly exceeds the other parts by a factor of 
~4(p/nmc 2 ) in the outflow channel. Since the enthalpy 
flux carries larger energy per cross section, the system 
can sustain faster energy input from the inflow region. 
In other words, the enthalpy flux is a key to sustaining 
fast reconnection. 

The guide field introduced interesting changes to the 
reconnection structure. Since no process can annihilate 
the out-of-plane magnetic field, the inflows transport the 
upstream guide field to the narrow outflow region, and 



therefore the guide field is compressed inside the out- 
flow channel. This significantly changes the energy com- 
position in the reconnection outflow: the outgoing en- 
ergy flow is rather dominated by the Poynting flux of the 
compressed guide field. Although there is no theoretical 
proof, our simulation results conserved (or slightly in- 
creased) the ratio of the B^-Poynting flux to the plasma 
energy a y . If we employ a principle of a y conservation, 
we can crudely estimate that the ratio of the upstream 
B^-Poynting flux to the downstream By-Poynting flux 
and the downstream enthalpy flow is given by 1 : a Vi i n . 

The boundary layers around the central reconnec- 
tion jet exhibit a significant charge separation ("f p n p — 
Jen e ) / ("fpUp + j e n e ) ~ 0.5 to retain the electric field sys- 
tem. The relative motion of species is significant in these 
layers, too. We think these conditions are beyond the 
scope of the single-fluid RMHD approach. The present 
resistive RMHD codes solve the charge distribution p c 
and the current j, independently from plasma bulk mo- 
tion. Such an approximation is valid when the charge 
separation and the relative motion are small. In fact, 
in the boundary layers, we pointed out that the single- 
fluid RMHD pressure and the two-fluid pressure differ 
by a factor of ~ 1.4. We also note that Koide (2009) 
(see Section 7.1 for a summary) discussed the validity 
of generalized RMHD equations in detail. In our case 
of a pair plasma with relativistic pressure, the pressure 
condition (p p /p p ~ p e /p e ) is weakly violated and the 
"proper charge neutrality" (n p « n e ) breaks down in the 
boundary layers (Figure Bp) . 

We demonstrated that the fluid inertial effects are im- 
portant to sustain magnetic reconnection in the recon- 
nection region. This is a characteristic feature of the two- 
fluid model. By definition, a single-fluid RMHD model 
only takes care of the fluid bulk inertia, while it does not 
consider the inertial resistivity. The inertia effects violate 
the ideal MHD condition around the bifurcated bound- 
ary layers, too. In other words, inertial effects enhance 
an effective resistivity both in the reconnection region 
and near the discontinuities, playing a similar role as an 
anomalous-type resistivity. To mimic such physics, the 
nonrelativistic MHD simulations often employ resistiv- 
ity profiles based on the electric current, for example, 
V = Vo + JJi[max(0, j/p — v cr it)} 2 , but its relativistic ex- 
tension may not be straightforward. 
Using a single-fluid RMHD theory, 



Lyubarsky (2005) 



predicted that relativistic guide field reconnection in- 
volves a three-layer structure in the outflow region, with 
rotational discontinuities (outside) and slow shocks (in- 
side). Furthermore, his balance analysis yielded that the 
outflow flux is dominated by the compressed guide field 
flux between the rotational discontinuities. In this work, 
we showed a different three-layer structure in the out- 
flow region: the central hot outflow inside two charge- 
separated boundary layers. In our case, the rotational 
discontinuities cannot appear in the magnetically dom- 
inated upstream regions, because the system runs out 
of plasma to sustain the electric current. The upstream 
region does not have sufficient plasmas to sustain the 
charge separation for the vertical electric field E z , and 
therefore the B y -Poynting flux (—cE z B y /Air) is confined 
around the central dense channel. On the other hand, the 
other important prediction qualitatively holds true — the 
outward energy is dominated by the guide field Poynting 
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flux. Energy transfer from the Poynting energy to the 
plasma energy will take place somewhere in the down- 
stream edge, such as the interaction between plasmoids 
and the ambient medium. 

In summary, we carried out full two-fluid simulations 
of relativistic magnetic reconnection with a guide field 
and investigated the characteristic properties, such as the 
charge separation and the guide field compression. We 
demonstrated that the guide field drastically changes the 
composition of the output energy flux, from enthalpy- 
dominated flow to Poynting-dominated flow, potentially 
controlled by a y ^ n . Inertial effects play a role of the ef- 
fective resistivity, and so they violate an ideal frozen-in 
condition. Most importantly, we showed that the multi- 



fluid approach is very useful to study important rela- 
tivistic plasma problems. It will be important to further 
study the consistency between the single-fluid RMHD 
model and the multi-fluid model. 
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